<*>  . 


E.  Royce  Fletcher 
L  Gerald  Bowen 


r-~ 

\J  L.-'! 

fn  - 11  tU  fl {. 

\  U  -  it  tM  -i  &  HHH  i  i  i  i 


Technical  Progress  Report 

■  •  .  y^;-'  '  on  •  ■  ,V  .y 

Contract  No.  DA-41  -146-XZ-372 


D 


1 


Thha  wc^  ah  aspect  of  investigations  dealing  with 
f'4  tl&;&il|^icalTEiiects  of  Blast  from  Somhfe,  was 
'•  imjjpbrt^jiby  tfte  Defense"  Atomic  Support  Agency  of 
:;y  f|  the  Department  of  Defense. 

••.*  ^  ; y  y  "-  v  y. 

S&dn.in  whole  oi* in  part  is  permitted  for  ss# 

t ) 


.  r.'p . 


jAlbuquerqtief  New  Mexico 


rch 


J*  \~ 


tmmmm  m»  &• 


fostrifchtiaa  of  Shit?  haaxmtf:  S#:tfes25aS.i. 


.  iU> 


DASA  1859 


-•  iiiniiArn 

BLA&HMJUbtU 


yrjtMPi  &Tlftil£l  rfXfPK 
inAn^LH liunnL  uuu'w 


E.  Royce  Fletcher 
I.  Gerald  Bowen 


Technical  Progress  Report 
on 

Contract  No.  DA-49-146-XZ-372 


This  work,  an  aspect  of  investigations  dealing  with 
the  Biological  Effects  of  Blast  from  Bombs,  was 
supported  by  the  Defense  Atomic  Support  Agency  of 
the  Department  of  Defense. 

(Reproduction  in  whole  or  in  part  is  permitted  for  any 
purpose  of  the  United  States  Government. ) 

Lovelace  Foundation  for  Medical  Education  and  Research 
Albuquerque,  New  Mexico 


November  1966 


/ 


FOREWORD 


This  report  describes  experimental  and  theoretical  studies  of 
blast-induced  translational  effects  with  regard  to  dummies,  animals, 
stones,  concrete  building  blocks,  and  spheres.  Some  results  are  presented 
of  mathematical  model  studies  of  the  accelerative  as  well  as  the  decelerative 
phase  of  displacement.  The  present  study  represents  a  further  analysis 
oi  previously  reported  data  and  is  intended  to  clarify  some  of  the  problems 
associated  with  translation.  The  biological  implications  of  translational 
events  are  not  discussed  in  this  paper. 

This  report  and  some  additional  material  were  summarized  on 
October  6,  1966,  before  a  symposium  on  the  Biological  Effects  of  Air 
Blast  at  the  National  Academy  of  Sciences  sponsored  by  the  NAS-NRC 
Committee  on  Hearing,  Bioacoustics,  and  Biomechanics ,  "CHABA.  " 

On  October  13,  1966,  the  report  was  presented  at  the  New  York  Academy 
of  Sciences  conference  on  Prevention  of  and  Protection  against  Accidental 
Explosion  of  Munitions.  The  text  was  also  submitted  for  inclusion  in  the 
Annals  of  the  New  York  Academy  of  Sciences. 


ABSTRACT 


A  mathematical  model  was  developed  to  predict  the  time  displace¬ 
ment  histories  of  objects  translated  by  the  blast  winds  from  conventional 
or  nuclear  explosions;  these  predictions  were  then  compared  to  actual 
experiments.  The  objects  studied  varied  in  size  from  139  mg  spheres 
to  man  and  were  all  assumed  to  be  free. to  move  over  a  smooth  horizontal 
surface.  The  effects  of  ground  friction'  could  cither  be  included  or 
neglected,  but  when  they  were  considered  the  ground  friction  eventually- 
brought  the  objects  to  rest  after  the  winds  had  passed.  The  values  of 
ground  friction  used  were  determined  experimentally  and  were  found 
to  be  functions  of  the  velocity  and  mass  of  the  object  being  displaced. 

The  translational  model  was  general  enough  for  cither  classical 
or  non-classical  blast  waves  to  be  considered.  Results  for  a  chemical 
explosion  were  obtained  by  using  both  the  computed  blast  waves  of  various 
authors  and  the  experimentally  determined  blast  waves.  These  predicted 
results  were  compared  with  each  other  as  well  as  with  experimental 
data  obtained  with  steel  spheres.  The  model  was  used  to  determine 
dynamic-pressure  impulses  necessary  to  explain  the  measured  sphere, 
velocities  at  three  ranges  from  ground  zero. 

Another  mathematical  model  was  briefly  described  which  was 
developed  to  compute  the  detailed  two-dimensional  trajectories  of 
objects  as  they  roll,  slide,  and  bounce  along  the  ground.  The  model 
closely  predicted  the  measured  distances  between  bounces  and  the  total 
displacements  of  concrete  blocks  and  large  stones  and  thus  helped  to 
explain  the  mechanisms  of  tumbling  by  which  an  irregular  object  may 
become  airborne  during  both  the  accelerative  and  the  decelorative  phases 
of  displacement. 


*  Ground  friction 
dissipation  b^  sliding, 


is  used  here  to  designate  the 
rolling,  and  impacting  along 


process  of 
the  ground. 


energy 


ACKNOWLEDGEMENTS 


The  authors  wish  to  acknowledge  their  appreciation  as  follows: 
to  Clayton  S.  White  for  many  informative  discussions  and  for 
editorial  assistance  in  preparing  the  final  manuscript;  to  Robert  F.D. 
Perett  and  Walter  A.  Sclieurlc  for  their  assistance  in  gathering  and 
analyzing  the  data;  to  Mildred  G.  Elrick  and  Ray.  W.  Albright  for 
developing  some  of  the  computer  programs  used  in  this  report,  to 
Robeit  A.  Smith  and  Joyce  Blaine  for  preparing  all  the  illustrative 
material;  to  Mildred  E.  Blake  for  editorial  support;  to  Lee  Sanders 
and  I, aura  M.  Carrasco  for  editorial  and  secretarial  assistance;  and 
finally  to  the  Defense  Atomic  Support  Agency  of  the  Department  of 
Defense  for  contract  support. 


ill 


BLAST-INDUCED  TRANSLATIONAL  EFFECTS 


E.  Royce  Fletcher  and  I.  Gerald  Bowen 


INTRODUCTION 

The  Lovelace  Foundation  has  conducted  experiments  on  two 
operations  at  the  Nevada  Test  Site  (Bowen  ct  al.  ,  1956,  1963;  Taborelli 
et  al.  ,  1959)  and  on  two  operations  using  chemical  explosives  (Bowen 
et  al.  ,  1965;  Fletcher  et  al.  ,  1965a;  llowcn  and  Fletcher,  unpublished) 
in  order  to  evaluate  the  translational  effects  produced  by  the  blast 
waves  from  nuclear  and  conventional  explosions.  Experiments  were 
performed  on  five  nuclear  shots  ranging  in  yield  lrom  1 1  to  T3  kt  and 
on  two  high  explosive  (HE)  shots  having  yields  of  50  and  500  tons. 
Translational  velocities  and  in  some  cases  total  displacements  were 
determined  for  about  20,  000  objects,  such  .as  spheres,  animals, 
anthropomorphic  dummies,  stones,  concrete  building  blocks, 
window-glass  fragments,  and  steel  fragments. 

Theoretical  studies  and  laboratory  experiments  (Bowen  et  al.  , 

1961,  1962;  Fletcher  etal.  ,  1961;  Fletcher  et  al.  ,  unpublished)  were 
done  to  help  explain  the  mechanisms  involved  in  the  translation  of 
objects  by  blast  waves  and  thus  to  make  it  possible  to  predict  their 
behavior  when  exposed  to  any  specified  arbitrary  blast  wave.  The 
decelerative  phase  of  displacement  when  the  objects  come  to  rest  by 
sliding,  rolling,  and  bouncing  along  the  ground  has  also  been  studied 
both  theoretically  and  experimentally  (Fletcher,  1965b). 

In  the  present  paper  the  results  of  several  experiments  (mostly 
conducted  during  a  500-ton  TNT  shot)  involving  the  translation  of 
spheres,  animals,  dummies,  stones,  and  blocks  will  be  described. 
Theoretical  time -displacement  histories  for  the  translated  objects  were 
computed  using  a  revised  mathematical  model  of  blast-induced  translation 
(Fletcher  et  al.  ,  unpublished).  These  calculations  were  based  on  blast 
waves  predicted  by  various  authors  as  well  as  on  experimentally 
determined  blast  waves;  in  this  way,  the  various  predictions  can  be 
compared  not  only  to  experiment,  but  also  to  each  other. 

TRANSLATION  MODEL 

The  revised  mathematical  model  (Fletcher  et  al.  ,  unpublished) 
of  blast-induced  translational  effects  is  similar  to  the  earlier  models 
(Bowen  et  al.  ,  1961,  1962),  but  two  major  features  have  been  added; 

(1)  arbitrary  blast-wave  forms  can  be  used  and  (2)  the  effects  of  ground 
friction  can  be  included  if  the  translating  object  is  tumbling  over  smooth, 
flat  ground.  In  the  earlier  models  it  was  assumed  that  the  blast  waves 
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were  classical  in  form  and  that  the  objects  did  not  strike  anything  while 
they  were  being  translated.  Thus  the  total  displacement  of  a  translated 
object  could  not  be  computed  since  these  modcis  included  no  way  of 
describing  the  motion  of  the  object  during  decelcrative  tumbling 
along  the  ground. 

As  in  the  case  of  the  earlier  modcis,  time-displacement  histories 
are  calculated  with  the  aid  of  a  digital-computer  program  based  on  the 
mathematical  model.  This  program  will  accept  an  overpressure  wave 
and  either  a  dynamic -pressure  wave  or  a  wind-velocity  curve,  and  will 
generate  the  time -displacement  history  of  an  object  (whose  acceleration 
coefficient  may  vary  with  time)  which  may  either  be  assumed  to  be 
in  free  air  or  to  be  bouncing  and  sliding  along  the  ground  as  it  translates. 
If  an  equation  is  not  available  for  the  blast  wave,  a  "broken-line" 
approximation  may  be  used  and  in  this  way  precursor  and  other  non- 
classical  wave  forms  can  be  handled. 

The  effects  of  ground  friction  on  the  object  are  approximated  in 
the  model  by  an  average  decelcrative  force  which  is  assumed  to  act 
on  the  object  throughout  the  entire  motion.  The  value  of  this  average 
decelcrative  force  due  to  ground  friction  was  derived  from  empirical 
fits  to  data  obtained  by  dropping  animals,  stones,  and  concrete  blocks 
from  a  truck  traveling  at  various  speeds  and  measuring  the  total 
displacements.  The  force  per  unit  mass  (i.  e.  ,  the  deceleration)  was 
found  to  be  a  function  of  velocity  and  mass  for  animals  and  a  function 
of  velocity  only  for  stones  and  concrete  blocks.  Equations  for  the 
frictional  deceleration  will  b  •»  presented  in  the  next  section  along  with 
a  description  of  the  "calibration"  experiments  from  which  they  were 
derived. 

In  developing  the  model,  a  differential  equation  was  derived  which 
describes  the  motion  of  an  object  acted  upon  by  blast  winds  and  ground 
friction.  Since  this  equation  could  not  be  solved  explicitly,  a  linear 
approximation  to  the  solution  over  finite  time  steps  small  enough  to 
introduce  negligible  errors  in  the  final  answers  was  used  in  the  computer 
program.  The  computed  acceleration,  velocity,  and  displacement  of  the 
object,  as  well  as  the  pertinent  blast-wave  parameters,  are  plotted  and 
printed  out  at  any  specified  interval  of  time. 

CALIBRATION  EXPERIMENTS 

In  order  to  use  the  translation  model  to  predict  the  complete  time- 
displacement  histories  of  objects  bouncing  along  the  ground;  it  was 
necessary  to  include  the  effects  of  ground  friction  in  the  model  in  some 
way.  Therefore,  experiments  to  study  ground  friction  (Fletcher,  1965b) 
were  performed  in  which  more  than  500  objects  were  dropped  on  smooth 
flat  ground  (a  graded  airstrip  on  an  alluvial  plain)  from  a  truck  traveling 
at  speeds  ranging  from  10  to  60  mph.  Total  displacements  of  the  objects 
were  recorded  as  well  as  the  distances  between  individual  impacts  if  the 
imprints  could  be  located.  The  experimental  objects  included  stones 


*Defined  as  the  area  presented  to  the  wind  multiplied  by  the  drag 
coefficient  and  divided  by  the  mass. 
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(1.3  -  116  lb),  ordinary  concrete  building  blocks  (26,  0  -  32.  5  lb), 
concrete-filled  building  blocks  (55.0  -  56.  5  lb),  guinea  pigs  (1.2  -  1.7  lb), 
rabbits  (5.2  -8.1  lb),  dogs  (38.0  -  44.6  ib),  and  a  goal  (59  lb).  The 
concrete  building  blocks  (7.  5  x  7.  5  x  15.  5  inches)  and  the  animals, 
which  had  been  sacrificed  by  injections  of  Nembutal,  were  dropped  in 
various  orientations  to  determine  variations  in  the  subsequent  motion. 

* 

The  animal  displacements  showed  no  variation  with  orientation 
but  did  show  a  significant  variation  with  mass;  for  a  given  speed  the 
larger  animals  required  a  slightly  longer  distance  to  stop  than  did  the 
smaller  ones.  A  dimensional  analysis  approach  was  used  in  an  attempt 
to  scale  from  one  animal  to  another  and  thus  to  explain  the  variation  in 
total  displacement  with  mass.  It  was  assumed  that  the  total  displacement 
of  mimnls,  S,  was  only  a  function  of;  (1)  the  initial  velocity,  2.  0  ; 

(2)  the  acceleration  of  gravity,  g;  (3)  the  mass  of  the  animal,  m; 

(4)  the  average  density  of  the  animal,  p;  (5)  the  coefficient  of  sliding 
friction,  (6)  the  coefficient  of  rolling  friction,  pr.  Then, by 

applying  the  principles  of  dimensional  analysis  and  assuming  that  g, 
p,  andp.s  are  constant  for  all  the  animal  drops,  the  following  relation 
was  derived; 

S  /  m  =  f(  h  /  m  ,  p  /  m  ) . 

It  should  be  noted  that  dimensional  analysis  cannot  predict  the  nature 
of  the  undetermined  function,  f.  Therefore,  in  an  attempt  to  determine 
f  empirically,  a  plot  was  made  ofS/m1-/^  vs  Zj  m1'  .  On  this  plot, 
the  datum  points  for  all  four  species  fell  along  a  single  curve  (a  separate 
curve  was  needed  for  each  species  on  a  plot  of  S  vs  Z^  indicating  that, 
within  experimental  accuracy,  S/m1'  ^  is  only  a  function  of  Z  /ni!'”and 
not  a  function  of  u^/m1'  .  This  could  be  interpreted  in  either  of  two 

ways:  (1)  p  r  /  m 1 '  ^  is  a  constant  or,  at  least  if  it  varies,  the  effect  on  S 
was  too  small  to  measure  in  the  calibration  experiments,  or  (2)  pr  should 
not  have  been  included  in  the  list  of  parameters  of  which  S  is  assumed  to 
be  a  function.  It  was  further  noted  that  when  the  plot  of  S/m1  ^  vs  ""Zj m1 '  ^ 
was  made  on  log-log  paper,  the  single  curve  appeared  to  be  a  straight, 
line. 

In  view  of  the  above  considerations,  it  was  decided  to  analyze  the 
animal  data  in  the  form 

log1Q  [S  (165/m)1/3j  =  A  +  B  log1Q  [ Zq(  1 65/m) 1 /6 ] 

where  A  and  B  are  constants  and  165/m  is  used  instead  of  1/m  so  that 
the  chart  prepared  to  show  this  relation,  FIGURE  1,  can  be  read  directly 
as  S  vs  ZQ  for  a  165-lb  animal.  For  an  animal  of  any  other  mass,  some 
calculation  would  be  involved  to  obtain  S  as  a  function  of  ZQ  from  the 
chart.  From  the  230  datum  points  in  FIGURE  1,  it  can  be  seen  that  the 
results  for  all  four  species  are  in  good  agreement  and  that  the  scatter  is 
quite  low,  the  standard  error  of  estimate  in  S/m1'  3  being  approximately 
1 1  percent. 

*  Regardless  of  the  original  orientation,  the  animal  assumed  a 
rolling  position  about  its  long  axis. 
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Since  tlic  equation  in  FIGURE  1  gives  the  scaled  average  stepping  distance 
as  a  function  of  the  scaled  initial  velocity  ,  it  is  possible  by  differentiating 
to  obtain  the  incremental  increase  in  stopping  distance  associated  with  an 
incremental  increase  in  initial  velocity  and  thereby  to  obtain  the  average 
deceleration  due  to  ground  friction.  JX  =  Z0  d(Z^)/dS,  as  a  function  of 
velocity  and  mass.  The  deceleration  Found  in  this  manner  was 


D  =  8.9029  [Z  (165/m)1'76]0- 38308 

•  7,  O  2 

where  Z0  is  i:i  ft/scc,  m  in  lb,  and  Da  in  ft/sec  .  From  this  deceleration 
it  is  possible  to  compute  the  time  required  for  the  animal  to  conic  to  rest. 
A  few  of  these  stopping  times  are  indicated  in  FIGURE  1  on  tile  scale 
labeled  t(165/)w^.  Thus,  the  slopping  times  can  be  read  directly 
for  a  165 -lb  animal  but  must  be  computed  for  an  animal  of  any  other  mass. 
In  order  to  facilitate  this  cirnpula ti on,  the  equation  for  t  as  a  function 

ol  Z  and  m  is  also  given  in  FIGURE  1. 
o 


Unlike  the  animal  displacements,  those  for  the  stones  showed  no 
variation  with  mass  over  the  range  of  the  stone  sizes  tested  (1.3  to  116  lb); 
i.e.,  stone  displacements  are  a  function  of  initial  velocity  but  not  of 
mass.  The  concrete -block  displacements,  on  the  other  hand,  did  vary 
with  orientation,  again  unlike  those  for  animals,  but  when  the  block 
data  for  all  orientations  were  combined,  there  was  no  significant  difference 
between  either  the  hollow-  or  the  filled-block  displacements  and  the 
stone  displacements.  Consequently,  in  order  to  predict  the  displacements 
of  cither  stones  or  blocks  dropped  in  random  orientations,  all  the  stone 
and  block  data  were  analyzed  together  as  shown  in  FIGURE  2.  The  datum 
points  represent  the  geometric  means  of  the  particular  group  of  blocks 
or  stones.  The  number  of  individual  experiments  in  each  group  is 
indicated  near  the  plotted  points.  (Individual  rather  than  grouped  data 
were  used  in  the  regression  analysis.  )  The  spread  in  the  data  is  so 
large  that  the  total  displacements  for  only  about  68  percent  of  similar 
individual  experiments  can  be  predicted  with.  41  percent.  The  stone 
and  block  displacements  have  a  greater  variability  than  the  animal 
displacements  because  the  motion  of  the  stones  and  blocks  is  not  as 
uniform  as  that  of  the  animals.  That  is,  animals  always  assume  the 
same  rolling  orientation  and  tend  to  stay  close  to  the  ground  whereas 
the  stones  and  blocks  are  mostly  airborne,  the  energy  dissipated  upon 
impact  being  greatly  dependent  on  orientation  as  well  as  velocity.  It 
should  be  noted  that  these  objects  were  dropped  from  a  truck  at  a  height 
of  approximately  2.6  ft  resulting  in  a  vertical  impact  velocity  of  about 
13  ft/sec  in  addition  to  their  horizontal  velocity',  Zc>.  The  deceleration 
for  stones  and  blocks,  Dg,  was  computed  (as  in  the  case  of  the  animals) 
and  found  to  be 


D  =  15.40  T Z  I0, 0996 

t»  O  2 

where,  as  before,  Z0  is  in  ft/scc  and  Dg  in  ft/sec  .  A  few  stopping 
times,  t,  and  the  formula  for  computing  them  are  also  indicated  in 
FIGURE  2  and  were  determined  in  the  same  manner  as  that  used  for 
computing  the  stopping  times  for  animals. 
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The  ai  imal  displacements  obtained  in  the  "calibration"  experiments 
have  beer,  compared  (Fletcher,  1965b)  to  the  results  obtained  by 
Anderson  et  al.  (1961)  where  goats  and  dummies  were  blown  from  an 
open-end  shock  tube  and  tumbled  to  a  stop  over  a  "grassy  pasture.  " 
Reasonably  good  agreement  was  found  between  the  results  of  the  two 
experiments  considering  the  variations  in  the  experimental  techniques 
and  the  inherent  differences  between  animals  and  dummies.  The 
experimental  stone  and  block  displacements  have  not  been  compared 
to  the  results  obtained  by  any  other  experimenters. 


SPHERE  TRANSLATION 


During  Operation  SnowBall  an  experiment  was  conducted  to  measure 
for  steel  spheres  of  different  sizes  their  blast-induced  velocities 
approximately  maximum  value  (Fletcher  et  al.  ,  1965a).  These  velocities 
could  then  be  related  by  means  of  the  acceleration  coefficient  to  those 
measured  photographically  for  goats  and  anthropomorphic  dummies  in 
another  experiment  to  be  described  later.  The  sphere  velocities  were 
also  used  to  estimate  the  dynamic-pressure  impulse  at  each  of  the  three 
ranges  where  they  were  placed. 


Impact  velocities  were  determined  by  trapping  the  spheres  in  two-inch 
layers  of  expanded  polystyrene  (Styrofoam)  and  later  measuring  their 
depths  of  penetration.  Two  pieces  (each  1  ft  x  3  ft  x  2  in)  of  Styrofoam 
were  used  at  each  station  placed  at  965-,  790-,  and  580-ft  ranges  where 
the  predicted  peak  overpressures  were  10,  15,  and  30  psi,  respectively. 

These  installations  were  identified  as  10SB,  15SB,  and  30SB,  the 
number  referring  to  the  predicted  overpressure  and  SB  representing 
SnowBall.  Each  sheet  of  Styrofoam  was  cemented  in  a  vertical  position 
to  a  piling  set  in  concrete  with  the  bottom  of  the  sheet  at  ground  level 
(FIGURE  3).  The  spheres  were  located  4.0  ft  upstream  from  the 
Styrofoam  at  heights  of  0.  8,  1.  8,  and  2.  8  ft  above  the  ground.  The 
initial  distance  between  the  spheres  and  the  Styrofoam  was  chosen  from 
theoretical  considerations  to  trap  the  spheres  at  approximately  their 
maximum  velocity.  The  spheres  were  placed  with  the  centers  two 
diameters  apart  in  shallow  holes  (0.  1  diameter  deep)  drilled  in  steel 
bars  which  were  one  inch  wide  and  one-fourth  inch  thick.  These  bars 
were  welded  to  vertical  steel  bars  (FIGURE  3)  set  in  the  ground. 

The  geometric  means  and  the  geometric  standard  deviations  of  the  -y 

impact  velocities,  V,  (as  computed  from  the  depths  of  penetration)  are 

shown  in  FIGURE  4  as  a  function  of  the  acceleration  coefficient,  a,  of  t 

the  spheres.  The  curves  are  the  predicted  V-vs-a  relations,  obtained 

using  the  blast  waves  from  the  indicated  references  and  the  mathematical 

model  of  blast-induced  translational  effects.  When  this  model  was  used 

to  compute  the  sphere-impact  velocities,  it  was  assumed  that  the  spheres 

were  in  free  air  (since  they  do  not  impact  the  ground  in  translating  from 

their  mount  to  the  Styrofoam)  and  that  their  acceleration  coefficients  were 

those  based  on  a  drag  coefficient  of  0.  47  (Hoerner,  1958).  The  calculations  ! 

were  discontinued  after  the  spheres  had  traveled  four  ft,  the  point  at  ’ 

which  they  should  impact  the  Styrofoam. 
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The  curves  based  on  Brodc's  work  (  1957)  were  obtained  by  finding 
the  reflection  factor  for  a  surface  burst  which  would  predict  the  measured 
peak-overpressure-vs -range  relation  from  SnowBall.  By  using  this 
reflection  factor  (2.  55)  along  with  the  ambient  pressure  (13.67  psi) 
and  ambient  speed  of  sound  (1142  ft/scc),  the  complete  set  of  blast- 
wave  parameters  could  be  obtained  as  functions  of  time.  The  curves 
computed  with  these  parameters  are  somewhat  high  for  15SB  and  30SB, 
but  seem  to. correspond  welbto  the  data  at  10SB.  Part  of  the  discrepancy 
may  be  explained  by  the  fact  that  an  unrealistically  high  reflection 
factor  of  2.  55  was  necessary  to  obtain  the  correct  overpressure  vs  range 
relation  whereas  a  reflection  factor  of  1. 63  was  calculated  for  SnowBall 
by  Keefer  (1965).  This  high  reflection  factor  (i.  e. ,  high  yield)  led  to 
longer  durations  and  correspondingly  greater  dynamic-pressure  impulses 
and  thus  predicted  velocities  which  are  too  great. 

The  curves  based  on  Dewey's  work  (1964)  fit  the  data  much  better, 
particularly  at  the  30SB  station,  although  the  predictions  at  15SB  and 
30SB  still  appear  to  be  somewhat  high.  The  better  fits  are  not  surprising 
in  view  of  the  fact  that  Dewey's  formulas  are  not  strictly  theoretical 
but  rather  are  empirical  fits  to  data  gathered  on  smaller  but  similar  shots 
(i.  e.  ,  hemispheres  of  TNT  placed  on  the  ground).  The  deviations  from 
the  experimental  velocities  may  result  from  the  fact  that  Dewey's  formulas 
are  not  realistic  for  times  beyond  the  limit  of  the  experimental  data 
since  they  predict  the  winds  will  remain  positive  forever;  i.  e.  ,  the  wind 
speed  reaches  a  minimum  value  and  then  starts  increasing  without  having 
gone  negative.  Although  some  of  our  experiments  (see  following  section) 
actually  did  show  evidence  of  an  "extended"  duration  at  the  10SB  station, 
no  such  evidence  was  seen  at  the  30SB  station  of  Operation  SnowBall. 

The  velocity  vs  a  curves  based  on  the  work  of  Mills  et  al.  (1964) 
are  high  at  all  three  stations.  However,  their  data  do  predict  the 
experimental  peak  overpressure-vs-range  relation  for  a  reflection  factor 
of  1. 63,  and  their  predicted  durations  for  the  positive  overpressure  are 
within  10  percent  of  those  measured  by  Keefer  (1965).  It  would  thus 
appear  that  the  dynamic -pressure  impulses  of  Mills  et  al.  are  too  large. 
Since  the  peak  dynamic  pressure  can  be  determined  from  the  peak 
overpressure  and  the  Rankine -Rugoniot  relations,  apparently  the 
dynamic -press ure  impulses  are  too  great  because  of  an  unrealistic 
shape  factor  or  duration,  or  both.  Note  that  the  predictions  based  on 
the  data  of  Mills  et  al.  and  Brode  are  in  agreement  for  the  30SB 
station  but  that  the  predictions  of  Mills  et  al.  become  progressively 
higher  than  those  of  Brode  as  the  overpressure  decreases. 

The  curves  labeled  "C"  in  FIGURE  4  were  computed  using  a 
combination  of  experimental  blast-wave  parameters  and  the  theoretical 
calculations  of  Brode:  i.  e.  .  the  exoerimental  peak  overpressure  and 
the  experimental  duration  of  the  positive  overpressure  waves  were  taken 
from  Minor  (1964)  and  Keefer  (1965),  the  peak  dynamic  pressure  was 
computed  from  theory,  the  dynamic -pressure  duration  was  calculated 
from  a  curve  of  (dynamic -pressure  duration)  /  (positive-overpressure 


duration)  vs  {peak  overpressure)  taken  from  Erode  (1957),  and  the 
overpressure  and  dynamic-pressure  impulses  were  determined  from 
curves  of  impulse  vs  overpressure  also  taken  from  Erode  and  based 
on  a  reflection  factor  of  1. 63.  The  measured  overpressure  curves 
appeared  to  be  nearly  classical  in  shape,  and  from  Dewey's  measurements 
(1964)  on  previous  shots,  it  was  anticipated  that  at  least  the  leading 
edges  of  the  dynamic -pressure  curves  would  also  be  near  classical. 

For  lack  of  more  detailed  information  the  overpressure  curves  were 
approximated  by  Friedlander  equations  with  linear  terms  and  the 
dynamic-pressure  curves  were  approximated  by  Friedlander  equations 
with  square  terms  as  is  done  in  The  Effects  of  Nuclear  Weapons 
(Glasstone,  ed.  ,  1962),  but  a  multiplicative  factor  was  introduced  into 
the  exponent  in  all  cases  so  that  the  desired  impulses  could  be  obtained. 

The  negative  phases  of  all  these  curves  were  not  used,  but  this  had  a 
negligible  effect  on  the  computed  impact  velocities  since  (1)  most  of  the 
spheres  impacted  before  the  negative  phase  arrived  and  (2)  those  spheres 
which  did  "see"  the  negative  phase  still  spent  most  of  the  time  traversing 
the  four -ft  distance  in  the  positive  phase  and  the  pressures  during  the 
positive  phase  were  much  larger  than  the  pressures  during  the  negative 
phase.  For  all  but  the  "C"  curves,  however,  the  negative  phases  were 
included  directly  as  they  were  computed  from  the  various  sources 
(note  that  Dewey's  dynamic -pressure  curves  did  not  have  a  negative 
phase  but  that  his  overpressure  curves  did). 

The  "C"  curves  seem  to  fit  the  data  reasonably  well  indicating  that 
Brode's  impulses  based  on  a  realistic  reflection  factor  are  reasonable 
even  though  the  corresponding  predicted  overpressure  vs  range  relation 
is  not.  The  computed  velocities  are  much  more  dependent  on  dynamic - 
pressure  impulse  than  on  duration;  hence,  the  good  agreement  between 
the  predicted  and  experimental  sphere  velocities  does  not  necessarily 
mean  that  the  predicted  dynamic-pressure  durations  are  close  to  the 
experimental  values. 

In  a  previous  paper,  (Fletcher  et  al.  ,  1965a)  dynamic -pressure 
impulses  were  estimated  by  the  following  procedure:  The  impact  velocity 
vs  a  data  at  each  station  were  approximated  using  a  least-squares 
regression  analysis  of  the  form 

V  =  Kj  o  t  K2  cZ 

where  Kj  and  are  constants.  It  was  shown  that  the  dynamic-pressure 
impulse  can  be  taken  to  be  equal  to  the  regression  coefficient  Kj.  The 
impulses  estimated  in  this  way  appear  in  Table  I  under  the  column 
labeled  "estimated  from  spheres.  "  The  impulses  derived  from  Brode 
using  a  reflection  factor  of  1. 63  and  used  in  computing  the  "C"  curves 
in  FIGURE  4  are  also  listed  in  Table  I  under  the  column  labeled  "predicted 
by  Brode."  Note  that  these  two  estimates  of  the  impulses  are  in  fair 
agreement,  the  biggest  difference  being  at  the  10SB  station.  The 
dynamic-pressure  impulses  of  Brode  were  adjusted  (leaving  all  other 


11 


/ 


TABLE-  I 


DYNAMIC-PRESSURE  IMPULSES  IN  PSI-MSEC 


Station 

Estimated- 
from  Spheres 

<= 

Predicted 
by  Erode  ( 1957) 

Adjusted 

Impulses 

10SB 

140 

120 

130 

1 5SB 

170 

190 

170 

30SB 

340 

340 

340 

*  Using  a  reflection  factor  of  1. 63.  These  impulses  were  used 
in  deriving  the  "C"  curves  in  FIGURE  4. 

Adjusting  the  impulses  to  obtain  curves  which  fit  the  data  in 
FIGURE  4  better  than  the  "C"  curves. 


quantities  including  peak  values  and  durations  the  same  as  used  in 
computing  the  "C"  curves)  to  find  the  impulses  which  best  fit  the  data 
plotted  in  FIGURE  4.  These  "adjusted"  impulses  are  also  listed  in  Table 
1  under  the  last  column.  Note  that  for  155B  and  30SB  the  values  agree 
with  the  estimates  from  the  regression  analysis,  but  for  10SB  the  adjusted 
value  is  somewhat  lower  than  the  impulse  estimated  from  the  sphere 
data.  This  discrepancy  apparently  results  from  the  fact  that  at  the 
10SB  station  the  mean  velocity  of  the  spheres  with  an  a  of  0.  139  fm/lb 
is  lower  than  what  would  be  predicted  based  on  the  velocities  of  the  other 
spheres  captured  at  this  station.  This  lower  velocity  would  tend  to 
increase  the  impulse  estimated  by  the  regression  analysis  and  decrease 
the  impulse  estimated  by  the  adjusted  Erode  analysis. 


The  close  agreement  of  these  impulses  and  the  moderate  scatter 
in  the  data  indicate  that  the  dynamic -pressure  impulses  can  be  approximated 
within  the  accuracy  of  the  data  by  using  the  simplified  regression  analysis 
rather  than  using  the  translation  model.  This  is  a  faster  and  simpler 
approach  which  still  allows  dynamic -pressure  impulses  from  short- 
duration  blast  waves  to  be  computed  to  a  reasonable  degree  of  accuracy 
without  the  details  of  the  blast-wave  parameters  being  considered  as 
functions  of  time. 


ANIMAL  AND  DUMMY  TRANSLATION 


During  Operation  SnowBall  an  experiment  (Bowen  et  al.  ,  1965) 
was  also  conducted  to  determine,  with  motion  pictures,  displacement 
and  velocity  as  functions  of  time  for  goats  and  anthropomorphic  dummies 
placed  at  the  same  ranges  as  the  10SB,  15SB,  and  30SB  sphere  stations. 


12 


_T*1  * 


- 


»«fiss<Stv 


Other  objectives  were  to  measure  the  total  displacements  of  the 
experimental  objects,  to  study  the  tumbling  mechanisms  during 
displacement,  and  to  compare  the  measured  velocities  with  those 
predicted  with  the  spheres. 

At  each  station  two  dummies  and  one  goat  were  placed  on  an 
asphalt  pad  on  which  a  grid  had  been  painted  and  which  was  lined 
on  cither  side  with  stripped  poles,  both  serving  as  reference  lines 
in  the  field  of  view  of  two  cameras  mounted  on  1  8 -ft  towers  placed  on 
cither  side  of  the  pad.  The  anthropomorphic  dummies  used  in  this 
experiment  weighed  165  lb  and  their  joints  were  adjusted  to  give 
sufficient  tension  for  the  dummy  to  stand  erect.  One  dummy  was  standing, 
facing  away  from  the  blast,  and  the  other  dummy'  was  prone,  face 
down,  with  its  head  toward  the  blast.  The  vertical  dummy'  was  supported 
by  leaning  it  approximately  2°  from  the  vertical  against  a  horizontal 
section  of  pipe  at  about  the  level  of  the  small  of  the  back.  The  goats, 
which  had  been  sacrificed  with  an  overdose  of  anesthesia  about  three 
hours  before  the  shot,  were  supported  in  a  normal  standing  position 
side-on  to  the  winds  with  the  aid  of  a  rack  designed  in  such  a  way  that 
translation  was  unhindered.  Thus,  both  the  goats  and  the  standing 
dummies  presented  a  maximum  area  to  the  oncoming  blast  wave. 

To  aid  in  the  interpretation  of  the  film  records,  the  positions 
at  various  times  of  the  initially -standing  dummy  and  the  goat  were 
sketched  on  a  single  drawing  of  the  asphalt  pad  as  shown  in  FIGURE  5 
for  the  10SB  station.  The  dummy  moved  out  feet  first  without  any 
apparent  rotation  about  its  long  axis.  For  a  time  the  resulting  airfoil 
tended  to  keep  the  dummy  airborne,  the  first  impact  being  on  the  dummy's 
head  after  about  547  msec  and  1 5 -ft  displacement  as  indicated  by  outline 
(h)  in  b  IGURE  5.  Note  that  the  goat's  head  and  feet  initially  moved 
ahead  of  the  trunk  resulting  in  an  orientation  which  would  cause  the  blast 
winds  to  deflect  the  motion  of  the  goat  to  the  right  and  upward.  The 
effect  due  to  the  head  motion  took  precedence  o  .  er  that  due  to  the  foot 
motion  causing  the  goat  to  impact  on  its  feet  with  the  head  in  a  forward 
position  after  a  displacement  of  about  11  ft  as  shown  by  outline  (f) 
in  FIGURE  5.  Note  that  both  the  dummy  and  the  goat  deflected  to  the 
lett.  Ihis  is  a  result  of  non-radial  winds  whicli  were  clearly  visible 
on  the  films .  1  he  motion  of  the  initially  prone  dummy  is  not  shown 

because  it  was  rapidly  obscured  by  dust. 

The  dummies  and  goat  at  the  15SB  station  were  obscured  by  dust 
after  about  200  msec  but  until  that  time  they  hebaved  in  a  manner  similar 
to  those  at  10SB.  One  difference,  however,  was  that  at  15SB  the  standing 
dummy's  arms  were  blown  forward  thus  reducing  the  total  drag  area. 

At  the  30SB  station  all  experimental  objects  were  obscured  by  dust 
immediately  upon  arrival  of  the  blast  wave  so  that  no  analysis  of  the 
film  records  was  possible. 

The  complete  time -displacement  histories  of  the  goats  and  dummies 
were  computed  with  the  translation  model  (including  the  ground  friction), 
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using  the  same  blast  waves  that  were  assumed  in  computing  the  "C" 
curves  in  FIGURE  4,  since  these  waves  liau  predicted  the  sphere 
velocities  fairly  accurately.  Because  most  of  the  acceleration  phase  is 
over  before  much  rotation  has  occurred,  the  acceleration  coefficients 
of  the  standing  dummy  and  goat  were  assumed  to  remain  constant  at  their 
initial  values  of  0.  052  and  0.  046  ft^/lb,  respectively,  during  the  entire 
translation  in  order  to  simplify  the  calculations.  This  approximation  would 
not  be  justified  in  cases  where  the  duration  of  the  blast  wave  is  long  enough 
to  allow  dummies  and  goats  to  rotate  appreciably  before  the  winds  have 
passed  (Tabor elli  ct  al.  ,  1959).  Results  of  the  calculations  for  the  standing 
dummies  and  goats  arc  plotted  in  FIGURE  6  in  the  form  of  velocity-vs- 
displaccment  curves,  the  points  and  straight  lines  representing  the 
experimental  data.  The  straight-line  segments  represent  measured 
velocities  which  appeared  to  be  constant  over  a  range  of  several  feet 
beyond  a  displacement  of  about  two  feet.  The  datum  points  indicated  by 
circles  arc  the  measured  total  displacements  and  the  points  indicated 
by  squares  arc  the  velocities  after  a  displacement  of  four  ft  as  computed 
from  the  polynominal  fits  to  the  sphere  data  (Fletcher  et  al.  ,  19f>5a) 
making  use  of  the  initial  c's  of  the  dummies  and  goats. 

At  the  30SB  station  the  predictions  about  both  the  dummy  and  the 
goat  agree  reasonably  well  with  the  measured  total  displacements  and 
the  velocities  estimated  from  the  sphere  data  (no  analysis  of  the  film 
records  was  possible  at  this  station  because  of  the  dust).  At  the  15SB 
station  the  computed  and  measured  velocities  of  the  dummy  and  goat 
are  in  close  agreement  as  are  the  computed  and  measured  total 
displacements  of  the  dummy,  but  the  measured  displacement  of  the 
goat  (48  ft)  is  considerably  greater  than  the  predicted  displacement  (2  3  ft). 
At  the  1C3B  station  all  of  the  computed  velocities  and  displacements 
are  significantly  lower  than  the  experimental  values.  It  is  thus  evident 
that  in  three  out  of  six  cases  there  is  good  agreement  between  the 
experimental  data  and  the  computed  t: me -displacement  histories,  and  in 
the  three  other  cases  the  agreement  is  poor.  It  was  also  noted  that  in 
the  three  cases  where  the  agreement  was  good,  the  average  angular 
deviation  from  the  radial  (measured  from  the  initial  to  the  final  position) 
was  2.  5°  with  two  objects  deviating  to  the  left  and  one  to  the  right. 

However,  in  the  three  cases  where  the  agreement  was  poor,  the  average 
angular  deviation  from  the  radial  was  14°,  all  three  trajectories  being 
to  the  left.  This  suggested  that  the  anomalous  velocity-vs-displacement 
histories  were  in  some  way  connected  with  the  anomalous  non-radial 
winds  which  had  been  observed  in  the  pictures  taken  at  the  10SB  station. 

Several  other  observations  provide  further  evidence  for  the  presence 
of  an  anomalous  wind  pattern  at  10SB.  The  velocity  of  the  initially- 
standing  dummy  increased  from  30.  4  to  33.  ?  ft/sec  about  0.  5  sec  after 
the  blast  wave  arrived  even  though  the  dummy  was  in  tumbling  contact 
with  the  ground  by  this  time.  This  seems  to  indicate  the  presence  of 
positive  winds  long  after  the  main  overpressure  wave  had  passed.  The 
trajectory  of  the  dummy  was  initially  about  9°  to  the  left  of  a  radial 
line  but  must  have  changed  during  the  latter  portion  of  translation 


*  The  measured  duration  of  the  positive  overpressure  pulse  was 
0.23  sec. 
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since  it  was  found  19*  to  the  left.  It  might  be  argued,  however,  that  the 
additional  deflection  could  have  been  caused  by  the  random  nature  ol 
the  decelcrativc  tumbling.  The  total  radial  displacement  of  the  dummy 
■was  about  60  ft  although  the  corresponding  dummy  at  15SB  (175  ft  nearer 
to  ground  zero)  with  about  the  same  maximum  velocity  was  displaced 
only  34  feet.  The  behavior  of  the  prone  dummy  (which  was  placed  on 
its  face  with  the  head  toward  ground  zero)  at.  10SB  indicated  the  presence 
of  non-radial  winds.  Although  the  average  radial  displacement  was  only 
1.0  ft,  the  dummy  had  been  rolled  over  to  a  back-down  position,  the 
head  and  feet  having  moved  3.  5  ft  and  6.  3  ft  to  the  left,  respectively. 

In  an  effort  to  understand  the  anomalous  wind  pattern,  FIGURE  7 
was  prepared  showing  the  10SB,  15SB,  and  30SB  dummy  stations  and 
indicating  the  initial  and  final  positions  of  the  goats  and  dummies  as  well 
as  the  predicted  final  positions.  Although  the  anomalous  winds  probably 
did  not  have  a  sharp  boundary,  apparently  they  were  fairly  localized 
in  that  the  goat  at  thel5SB  station  suffered  a  much  greater  angular 
deviation  than  did  the  standing  dummy  (10.6*  as  opposed  to  3.0°) 
which  was  initially  only  8  ft  from  the  goat.  In  view  of  this,  a  line  was 
drawn  on  FIGURE  7  to  indicate  the  approximate  boundary  of  the  anomalous 
winds.  The  line  was  drawn  at  an  angle  of  14°  with  respect  to  the  radial 
since  this  was  the  average  angular  deviation  of  the  three  objects  presumed 
to  be  in  the  anomalous  winds.  Note  that  the.  U.S.  blast  line  was  situated 
well  outside  the  apparent  boundary  of  the  non-radial  winds  and  hence 
no  blast-wave  instrumentation  was  available  in  the  region  of  the  anomalous 
winds.  There  is  no  indication  that  any  of  the  sphere  stations,  which  in 
all  three  cases  were  located  more  than  200  ft  circumferentially  from 
the  dummy  stations  (FIGURE  7),  were  affected  by  anomalous  winds; 
however,  non-radial  winds  would  be  difficult  to  detect  since  the  distance 
of  travel  of  the  spheres  was  only  four  feet.  Unfortunately,  the  above 
analysis  does  not  indicate  the  causes  of  the  anomalous  winds  nor  the 
frequency  of  occurrence. 

STONE  AND  BLOCK  TRANSLATION 

In  1957  during  Operation  Plumbbob,  experiments  were  conducted  on 
two  nuclear  shots  to  study  stone  and  block  translations  (Bowen  et  al.  ,  1963). 
One  of  these  experiments  will  be  described  since  the  interpretations  of  the 
results  are  applicable  to  the  problem  of  translational  effects  from 
conventional  explosions.  On  Shot  Galileo  (11  kt,  500-ft  HOB),  a 
concrete-block  wall  was  constructed  at  a  ground  range  of  2750  ft  to 
study  the  breakage  and  displacement  of  the  blocks  due  to  the  blast  wave. 

This  wall  (FIGURE  8,  preshot  view)  which  weighed  about  four  tons  and 
contained  236  concrete  blocks  (7.  5  x  7.  5  x  16  inches,  34  lb)  was  located 
on  level,  graded  ground.  FIGURE  9  is  a  view  of  the  wall  after  it  had 
been  struck  by  a  near-classical  blast  wave  with  a  peak  overpressure  of 
8.  7  psi  and  a  duration  of  0.  58  seconds.  Note  that  many  of  the  blocks 
remained  intact  and  that  some  of  them  were  apparently  stopped  by  the 
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Figure  7.  Layout  of  the  three  dummy  stations  showing  the  initial  and 
final  positions  of  the  goats  and  dummies  as  well  as  the  predicted  final 
positions.  The  approximate  boundary  of  the  anomalous  winds  is  also 

indicated 
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three  Styrofoam  missile  traps  placed  behind  the  wall  (FIGURES  8  and  9). 
Final  resting  positions  for  the  larger  wall  fragments  (whole,  half,  and 
joined  blocks)  are  plotted  in  FIGURE  10.  One  block  (not  plotted  but 
indicated  at  the  top  of  the  chart)  was  found  as  far  as  403  ft  downstream 
and  170  ft  to  the  left  of  the  center  of  the  wall.  Note  the  absence  of  blocks 
behind  the  trap  installations. 

In  order  to  estimate  the  velocities  of  the  wall  fragments  from  their 
total  displacements,  it  was  necessary  to  use  the  relation  for  initial 
velocity  vs  displacement  for  stones  and  blocks  as  determined  from  the 
calibration  experiments  described  in  a  previous  section.  Further,  it 
was  necessary  to  make  three  simplifying  assumptions:  (1)  The  wall 
fragments  achieved  their  pe..k  velocity  before  they  were  displaced  any 
significant  distance.  This  is  approximately  true  for  objects  with  low 
acceleration  coefficients  such  as  those  for  the  larger  wall  fragments, 
<0.03  ft^/lb.  In  effect,  this  assumption  allows  us  to  neglect  the  blast 
wave  and  to  simply  think  of  the  wall  fragments  as  having  started  off 
with  an  "initial  velocity.  "  (2)  These  "initial  velocities"  were  horizontal 
(no  vertical  component)  and  followed  a  log  normal  distribution,  sometimes 
called  a  geometric  distribution.  (3)  All  the  blocks  from  each  of  the 
eight  layers  of  blocks  first  impacted  the  ground  at  the  point  predicted 
by  the  geometric  mean  "initial  velocity,  "  the  height  of  the  layer  above 
ground,  and  the  acceleration  of  gravity.  This  is  an  approximation 
since  the  actual  first  impact  points  would  be  scattered  about  this  distance 
due  to  variations  in  the  "  initial  velocities.  " 

With  the  above  assumptions  it  is  possible  to  compute  the  distribution 
of  displacements  if  a  geometric  mean  initial  velocity  and  a  geometric 
standard  deviation  are  assumed.  The  best  values  of  these  parameters 
were  determined  by  trial  and  error  and  the  results  of  this  analysis  are 
shown  in  FIGURE  11  where  the  downwind  distance,  d^,  is  plotted  against 
the  percent  of  the  wall  fragments  traveling  less  than  the  indicated 
distance.  The  circles  on  this  plot  represent  the  experimental  cumulative 
distribution  of  the  155  wall  fragments  and  the  curve  represents  the 
theoretical  distribution  for  a  geometric  mean  velocity  of  32  ft/sec  and 
a  geometric  standard  deviation  of  1.  55.  Since  these  two  parameters 
affect  the  curve  in  different  ways  (changing  the  velocity  tends  to  translate 
the  curve  while  changing  the  standard  deviation  tends  to  tilt  the  curve), 
both  are  closely  determined  by  the  process  of  fitting  the  data,  a  change 
of  as  little  as  4  percent  in  either  the  geometric  mean  initial  velocity  or 
the  geometric  standard  deviation  producing  a  curve  which  clearly  does 
not  fit  the  data  as  well  as  the  curve  shown  in  FIGURE  11.  However, 
because  of  the  approximations  and  the  large  scatter  in  the  data,  this 
docs  not  mean  that  the  true  values  lie  within  4  percent  of  the  estimates, 
but  it  does  show  that  the  possible  errors  introduced  by  the  process  of 
fitting  the  data  are  relatively  small. 

In  another  experiment  at  the  same  range  as  the  block  wall,  steel 
spheres  were  trapped  in  Styrofoam  at  approximately  their  maximum 
velocity.  It  was  shown  that  the  sphere  velocities  agreed  well  with  those 
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computed  with  the  earlier  translation  model  (Bowen  et  al.  ,  1963). 

Using  an  average  initial  acceleration  coefficient  (0.  0263  ft^/lb  assuming 
a  drag  coefficient  of  1.07)  for  the  blocks  in  the  wall,  a  peak  velocity  of 
30  to  33  ft/sec  would  be  predicted  from  the  sphere  data.  This  is  in 
excellent  agreement  with  the  geometric  mean  initial  velocity  of  32  ft/sec 
computed  from  the  block  displacements,  indicating  that  velocities  of 
blocks  from  walls  can  be  reasonably  estimated  by  assuming  the  free-field 
blast  wave  and  using  the  initial  acceleration  coefficient  of  the  blocks  — 
at  least  in  the  lower  pressure  regions  and  for  similar  types  of  block 
walls  which  are  not  reinforced.  There  is,  however,  the  possibility 
that  the  true  geometric  mean  velocity  was  a  little  higher  than  32  ft/sec 
since  some  of  t  ie  displacements  would  probably  have  been  greater 
if  the  traps  had  not  been  located  behind  the  wall.  This  effect  is  probably 
small  because  at  least  some  of  the  blocks  striking  the  traps  were 
fragmented  by  the  impact  and  hence  were  not  included  in  the  ''larger" 
fragments  used  in  this  analysis. 

DECELERATIVE  TUMBLING  MODEL 

A  mathematical  model  of  decelerative  tumbling  was  developed 
(Fletcher,  1965b)  to  study  the  mechanisms  of  tumbling  of  the  stones, 
blocks,  and  animals  and  to  be  able  to  extrapolate  to  greater  initial 
velocities  than  those  used  in  the  experiments.  In  this  model,  a  two- 
dimensional  rigid  object  is  assumed  to  slide,  roll,  lift-off,  and  impact 
along  perfectly  horizontal  ground  where  the  interactions  between  the 
ground  and  the  object  are  through  the  coefficients  of  sliding  and  rolling 
friction  and  the  coefficient  of  restitution  (assumed  to  be  zero).  The 
shape  chosen  for  the  translated  object  was  a  "generalized  ellipse"  of 
the  form  (x/a)n  +  (y/b)n  =  1  where  1  <  n  <  co.  By  choosing  appropriate 
values  of  a,  b,  and  n,  the  object  can  be  made  to  approximate  rectangular 
blocks  or  rocks  and  animals  with  various  degrees  of  roundness. 

It  is  of  interest  to  note  the  mechanisms  of  the  model  by  which 
objects  lift  off  the  ground,  particularly  since  the  coefficient  of  restitution 
between  the  ground  and  the  object  is  assumed  to  be  zero  and  therefore 
the  objects  do  not  bounce  in  the  usual  sense  of  the  word.  When  the  force 
(frictional  and  normal)  of  the  ground  against  the  object  is  in  such  a 
dir  ection  as  to  result  in  an  upward  component  of  velocity,  the  object 
can  lift  off.  In  general  then,  without  the  sliding  frictional  force,  the  object 
would  slide  and  roll  (rather  than  bounce)  until  it  came  to  rest  because 
of  the  rolling  friction  losses. 

Because  of  the  complexity  of  the  mathematical  model,  computations 
were  performed  on  an  electronic  digital  computer.  FIGURE  12  shows 
some  of  the  computed  output  for  an  object  with  major  and  minor  axes 
of  1 . 0  and  .05  ft,  respectively,  (n  =  3)  which  had  an  initial  horizontal 
velocity  of  100  ft/sec,  an  initial  downward  velocity  of  10  ft/scc,  and  no 
initial  angular  velocity  (the  physical  constants  are  indicated  on  the 
figure).  The  height  of  the  center-of-mass,  angular  velocity,  distance 
traveled,  and  energy  per  unit  mass  are  all  plotted  as  functions  of  time. 


Figure  12.  Some  of  the  translation  parameters  computed  using 
the  decelerative  tumbling  model.  The  shape  of  the  object  being 
displaced  as  well  as  its  initial  velocity  and  orientation  are  indicated 
in  the  lower  fourth  of  the  figure. 
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Note  that  the  velocity  and  energy  are  "discontinuous"  when  an  impulsive 
impact  occurs.  When  the  object  is  airborne,  the  height-above-ground 
curve  should  be  parabolic  rather  than  triangular  as  plotted  by  the  computer 
for  the  sake  of  simplicity;  however,  the  peak  heights  as  shown  in  FIGURE 
12  are  correct.  Note  that  the  object  is  airborne  most  of  the  time,  this 
being  typical  of  the  actual  motion  of  stones  and  blocks  but  not  of  the 
rounder  softer  animals. 

The  first  detailed  calculations  made  with  the  mathematical  model 
were  designed  to  simulate  the  motion  of  the  concrete  blocks  dropped 
from  the  moving  truck  and  described  in  a  previous  section.  Values  of 
rolling  and  sliding  friction  coefficients  were  estimated  and  subsequently 
adjusted  in  order  to  achieve  reasonable  agreement  with  the  experimental 
results.  During  the  experiments  eight  blocks  were  dropped  at  29.  3  ft/sec 
with  their  long  axis  perpendicular  to  the  motion,  four  dropped  on  edge, 
and  four  flat  on  one  face.  FIGURE  13  shows  both  the  computed  and  the 
measured  distances  as  a  function  of  bounce  number.  To  introduce  scatter 
into  the  theoretical  data,  five  calculations  were  made  for  each  orientation; 
these  runs  assumed  the  block  was  dropped  at  exactly  the  desired  angle 
and  at  ±5*  and  ±10*  from  the  desired  angle.  Notice  that,  after  the  first 
bounce,  both  the  theoretical  and  experimental  dafa  divide  into  two  groups 
according  to  orientation.  The  curve  in  FIGURE  13  separates  both  the 
theoretical  and  experimental  distances  associated  with  one  orientation 
from  the  same  values  associated  with  the  other  orientation  (after  the 
first  bounce) . 

FIGURE  14  is  a  plot  of  initial  horizontal  velocity  vs  total  distance 
traveled  for  all  the  blocks  dropped  in  the  two  orientations  shown  in  FIGURE 
13.  The  53  experimental  drops  are  approximately  divided  evenly  among 
six  velocities:  14.  7,  29.  3,  44.  0,  58.  7.  73.  3,  and  88.  0  ft/sec.  Again 
the  scatter  in  the  theoretical  runs  was  obtained  by  assuming  the  blocks 
were  dropped  at  exactly  the  desired  angle  and  at  angles  deviating  from 
this  angle  by  ±5*  and  ±10?  Thus  five  runs  were  made  for  each  orientation, 
or  a  total  of  ten  theoretical  runs  for  each  velocity.  Similar  to  those  in 
FIGURE  13,  the  theoretical  points  in  FIGURE  14  are  slightly  displaced 
on  the  ordinate  so  they  can  be  easily  distinguished  from  the  experimental 
points.  Note  that  neither  the  theoretical  nor  the  experimental  points 
are  separated  by  orientation  as  they  were  in  FIGURE  13. 

From  FIGURE  13  and  FIGURE  14  it  can  be  seen  that  the  theoretical 
and  experimental  results  for  blocks  agree  closely  in  average  displacement 
and  scatter  in  the  displacements.  Stone  displacements  can  be  computed 
by  averaging  the  total  displacements  associated  with  all  orientations  of 
the  blocks  as  is  indicated  in  FIGURE  2.  At  this  writing,  no  attempt  has 

hf>f*n  m ar!p  tn  rnmnnfp  t Vi Hprolpratiuo  T'rn/,fl£S  f^r  animals. 
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DISCUSSION 

In  addition  to  the  translational  effects  discussed  in  this  paper,  there 
are  several  other  related  areas  of  interest  which  are  important,  but  have 
not  been  included  because  of  space  limitations.  Some  of  these  effects 
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should  at  least  be  mentioned  even  though  it  is  beyond  the  scope  of  this 
paper  to  discuss  the  damage  to  biological  or  physical  targets  that  might 
result  from  these  translational  phenomena. 

Of  particular  significance  is  the  problem  of  secondary  missiles, 
or  blast-energized  debris.  These  missiles  may  or  may  not.  represent 
a  biological  hazard  depending,  among  other  things,  on  the  velocities, 
sizes,  shapes,  and  materials  involved  If  the  blast-wave  parameters 
are  known  and  if  the  missiles  arc  relatively  free  to  move,  the  translation 
model  can  be  used  to  computc-the  time -displacement  histories  of  the 
objects  assuming  they  are  in  free  air  or  bouncing  along  the  ground. 
Additional  assumptions  are  sometimes  necessary  depending  on  the 
circumstances  of  exposure.  Appropriate  approximations  for  predicting 
the  time-displacement  histories  of  fragments  from  concrete -block  walls 
were  discussed  in  this  paper,  and  those  for  computing  the  motion  of 
window -glass  fragments  were  studied  by  Bowen  et  al ,  ($963). 

One  limitation  of  the  present  translation  model  is  the  assumption 
of  unidirectional  winds  which  produce  only  one -dimensional  displacements. 
It  has  been  shown  that  the  model  is  reasonably  adequate  for  computing 
displacement  associated  with  near-  classical  blast  waves  in  "clean" 
geometries  (i.  e.  ,  unperturbed  radial  winds),  and  although  verification  is 
not  now  available,  it  is  presumed  that  model  solutions  for  disturbed 
wave  forms  (precursor,  for  example)  would  be  satisfactory  provided 
appropriate  parameters  of  the  wave  could  be  specified. 

Other  problems  such  as  translation  in  shelters,  in  and  behind 
buildings,  and  in  forests  ,  do  not  fall  in  the  categories  mentioned  above. 

To  study  these  thoroughly,  it  would  be  necessary  to  develop  a  three- 
dimensional  translation  model  and  a  three-dimensional  dcceleiative 
tumbling  model.  Although  some  translation  experiments  have  been 
conducted  in  shelters  (White  ct  al.  ,  1957;  Richmond  et.  al.  ,  1959;  Bowen 
et  al.  ,  1963)  and  the  problems  have  been  approached  theoretically 
(Fletcher,  1965c),  much  more  experimental  and  theoretical  effort  is 
needed  in  this  field. 


f 


J; 

& 

£ 


! 


It  has  been  shown  that  the  translation  of  gravel  placed  in  front  of 
a  wall  at  a  distance  greater  than  the  structure's  height  or  width  was  not 
greatly  affected  by  the  presence  of  the  obstruction  (Bowen  et  al.  ,  1963). 

However,  an  adequate  theory  has  not  been  developed  for  the  translation 
of  objects  located  relatively  near  a  "large"  obstruction  when  the  blast 
wave  is  finite  in  duration,  although  a  theory  fo.i  fiat-top  blast  waves  has 

been  reported  (White  et  al.  ,  1965,  Fletcher,  1966c).  1 

1 

Another  translational  problem  associated  with  surface  bursts  involves 
crater  debris  or  ejecta.  Because  both  earth  and  air  motion  accelerate 
the  ejecta,  the  translation  model  cannot  be  used  to  predict  velocities 
and  displacements.  However,  ejecta  velocities  have  been  estimated  from 
the  distances  to  impact  with,  and  without,  the  arrival  times  (Bowen  et  al.,  1965).  } 
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SUMMARY 


A  revised  mathematical  model  of  blast-induced  translational  effects 
was  discussed  which  has  two  advantr »es  over  the  earlier  models:  (1) 
arbitrary  blast-wave  forms  can  be  vised  and  (d)  the  effects  of  ground 
friction  can  be  included  if  the  translating  object  is  tumbling  over 
smooth  flat  ground.  For  relatively  short-duration  blast  waves  (such  as 
those  produced  by  HE)  the  single  most  important  bla&l-wavo  parameter 
needed  to  predict  translation  is  the  dynamic -pressure  impulse. 

Experiments  involving  the  translational  velocities  of  steel  spheres 
were  described,  the  results  of  which  were,  used  to  estimate  dynamic- 
pressure  impulse  by  two  methods;  the  first  was  accomplished  by  means 
cf  the  mathematical  model  and  the  second  simplified  technique  required 
only  a  least-squares  analysis  of  the  experimental  data.  The  experimental 
sphere  velocities  wci  c  compared  to  the  velocities  computed  with  the 
translation  model  using  the  blast  waves  predicted  by  various  authors,  and 
the  discrepancies  were  discussed. 

In  order  to  understand  the  process  of  decelerative  tumbling,  a 
series  of  "calibration"  experiments  were  performed  in  which  concrete 
blocks,  stones,  and  animals  were  dropped  from  a  truck  traveling  at 
various  speeds  and  the  distances  between  bounces  and  the  total  displacements 
were  measured.  From  these  experiments  an  "average  ground  friction" 
term  was  derived  and  included  in  the  translation  model. 

The  model  was  tested  using  data  obtained  in  experiments  where 
dummies  and  goats  were  displaced  by  the  winds  from  a  000-ton  TNT 
detonation.  The  computed  velocities  and  displacements  were  in 
reasonable  agreement  with  the  experimental  results  except  in  those 
cases  where  anomalous  non-radial  winds  were  present.  The  model  could 
not  be  adequately  tested  for  the  latter  cases  since  appropriate  blast-wave 
parameters  were  not  available. 

An  another  experiment,  the  velocity  distribution  of  the  larger 
fragments  from  a  concrete -block  wall  blown  apart  by  the  winds  from  an 
11 -kt  nuclear  shot  was  reconstructed  from  the  measured  displacements 
of  the  fragments  and  the  "calibi  ntioa"  experiments.  The  geometric  mean 
velocity  of  the  fragments  was  shown  to  be  in  good  agreement  with  the 
predicted  velocity  for  a  single  loose  block  situated  in  fr  ;e  air. 

A  brief  description  was  given  of  a  mathematical  model  which  was 
developed  to  study  the  mechanisms  of  decelerative  tumbling  for  stones, 
blocks,  and  animals  and  to  be  able  to  extrapolate  to  greater  velocities 
than  have  heretofore  been  obtained  experimentally.  The  predictions  of 
the  model  were  shown  to  be  in  good  agreement  with  the  results  of  the 
"calibration"  experiments  for  block  translations. 
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